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1 Introduction 

Einstein field equations are nonlinear, and therefore cannot in general be solved 
exactly. By imposing the symmetry requirements of time independece and spa- 
tial isotropy, the Schwarzschild solution can be obtained, but this solution can 
not actually describe in the solar system because it is not static and isotropic. The 
Post-Newtonian approximation is a systematic method in which we can describe 
a system of slowly moving particles bound together by gravitational forces. In 
this method we will use an expansion of geometric objects as the metric tensor, 
in inverse powers of the speed of light using the parameter /3 = ^ . 

Several works had calculated the perihelion precession of planetary orbits based 
on Einstein's equations with the line element of Schwarzschild and the use of 
Euler-Lagrange equations. In this paper we present a simplified derivation using 
a post-newtonian approximation for a general spherically symmetric line ele- 
ment. In order to obtain the precession [T], we compare the Keplerian orbit with 
the curved space trajectory and making use of the invariance of Kepler's second 
law. 

2 Post-Newtonian Metric 

The Minkowski metric in polar coordinates is given by 



1310-0157 © 2005 Heron Press Ltd. 



1 



A. Larranaga, L. Cabarique 



ds^ = df - dr^ - r^d§^ - sin^ d dip^ 



(1) 



whereas a general curved metric with spherical symmetry has a line element that 
can be written as 



ds^ 



gtt{r)dr + grr{r)dr^ - rd'd^ - r sin^ Mip^. 



(2) 



In the post-newtonian approximation the metric components gu and g^r are ex- 
panded up to the fourth order in ^ = - in the form 



grr = -1 + + ffV, + O . 



(3) 
(4) 



Thus, the relation between flat and curved spaces for the radial and time coordi- 
nates can be considered as the transformation 



and 



dt' = 



dr' = 



1 (2) 1 (4) 

2 5 tt + 2 tt 



1(2) 



2 2 



1(4) 
9 r 



dt 



dr, 



(5) 



(6) 



where the binomial approximation has been applied. 
3 The Advance of the Perihelion 

In order to obtain the advance of the perihelion we will consider two elliptical 
orbits. First, the Keplerian orbit in flat space with coordinates {t, r) for which 
the element of area is given by 



dA= [ 
Jo 

This equation gives the second law of Kepler, 



rdrd'd = — di?. 
2 



dA _ 

'dt ~ ~2~di' 



(7) 



(8) 



The second elliptical orbit that we will consider is obtained in the coordinates 
it' , r') of the curved space. Here, the element of area is 
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dA' = / rdr'a-d 
Jo 



dA' = 



dA' 



R 




i?2 



1(2) 1(4) 

rr 



drd^ 



1 



R 



(2) (4) 
Q rr 9 rr 



dr di9 



(9) 
(10) 

(11) 



where we have used (|6]l. Kepler's second law is now 



dA' i?2 



dt' 



(2) (4) 
9 rr+ 9 rr 



(12) 



which, using equation (|5]l, gives 



dA^_ie_ 
dt' ^ 2 



1 - 



1 



9 rr [n + 9 rr [r) 



R^Jo 

Using again the binomial approximation it gives 

dA' i?2 ^ di) 



dr 



1 (2) 1 (4) 

1+^ 9 tt iR} + ^ 9 tt (R) 
(13) 



(14) 



di} 



where we have defined 



1 



R^J, 



9 rr in + 9 rr (r) 



dr 



1 1^'' ^RA 

(15) 



Applying this increase to change dd to dd' in a single orbit gives the expression 



dd' = Ad' 



<i>{R)di}. 



(16) 



In order to perform the integration in the right hand side we will consider the 
Keplerian ellipse using R {■&) = where / is the latus rectum and e is the 

eccentricity. 
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4 Particular Cases 

4.1 Schwarzschild Metric 

Up to second order in /3^, Schwarzschild's metric is 



ds^ = {l- —) dt^-(l + —] dr^ - r^d^^ - sin^ (17) 



r ) 



and therefore 



r J 



(2) , . (2) , ^ 2M 

9 rr = 9tt{n = — 



5 rr ('') = 9tt (r) = 0. 



Function $ (i?) is given by 



= 1 



r (r) dr 



1 



R J \ R J 
2,M 2M2 



R i?2 
and therefore the increasing Ai?' is 



Jo 



Ai?=27r r 



3M 



(18) 
(19) 



(20) 

(21) 
(22) 
(23) 



(24) 



In order to perform the integration, we will the Keplerian elUpse using R (i?) = 

TT-^ a. SO 



M'= [ 
Jo 



1 + ^ (1 + e COS 1?) + (1 + e cos 



2M2 
12 



di? (25) 



Ai9' = 27r + — — + — — (2 + e^) . 



(26) 
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Hence, the obtained perihelion advance for Schwarzschild's metric has the stan- 
dard value ^^j^ plus an additional term of order P^. 

4.2 Reissner-Nordstrm Metric 

Up to fourth order in 13, the metric for an electrically charged and spherically 
synmietric object is 



ds'= 1- 



2M 



dt^- 1 + 



2M Q2 



dr^-r'^d'd'^-r'^ sin^ Mif"^, 
(27) 



and therefore we will take the approximation terms 



(2) , . (2) 2M 

9 rr in = 9tt{n = — — 



(28) 



(4) , , 

9tt (r) = ^• 



(29) 



This time function $ (i?) is 



= 1 + 



R'' 

2M\ 



M 



= 1 + 



3M 2M2 g2 MQ 



R 



1 - 




M 




'r ~ 


2R'^ ) 










i?3 





1(4) 

2^* 



and the increasing Ai?' is now 



M' = 



'-L 



Ai>=27r r 



1 



3M 2M2 Q2 MQ^ 
~R^^^?~~2R^~ R3 



d'd. 



(31) 
(32) 
(33) 

(34) 



In order to perform the integration we consider again R [p) = xqij^j^. and 
therefore 



M' = 



'-1 



Ai?=27r 



1 + ^ (1 + ecosi?) 



"l2-"^j ^^ + "°^^) --^(1 + ecos^?) 

(35) 

5 



d-d 



A. Larranaga, L. Cabarique 



A^' = 27r+-^ + -^(2 + e2)_^(2 + ,2)__J^(2 + 3e2) 

(36) 

We obtained again the standard value ^^j^ and the order term for 
Schwarzschild's solution given above, but now we also obtain an additional term 
of order /3'^ solely due to the electric charge of the central body which agree with 
the perihelion precession reported in ||2l and |l3l. Finally we obtain a term that 
is proportional to the product MQ^ that is not reported in those papers. 

Finally, it is also important to note that the terms that contain the electric charge 
have a negative contribution to the perihelion precession. 

4.3 Schwarzschild-de Sitter and Schwarzschild-Anti de Sitter Metric 
The Schwarzschild-de Sitter metric is 



(37) 

where A is the cosmological constant and therefore we will take 

9 rr = 9 ttV') = — (38) 



^9rr {r) = ^Ou {r) = Ar^. (39) 
This time function <!> (i?) is 
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3M 2AP 5AMR 3AR^ 

~R ^ ~R2 4 — + 8 



M. (44) 



In order to perform the integration we consider again R{d) = y+tI^, and 
therefore 



A)?=2ir r 



3Af 

1 + — (1 + ecosi?) 



-^(l + ecos?9)^ 



SAAfZ 



2 74 



4(l + ecosi?) 4(l + ecosi?)^ 8 (1 + ecosi?)'' 
which for small eccentricities can be integrated as 



Ai?' = 27r + + (2 - 



I 



P 



birAMl SttAP 7rA2;4 



(46) 



where the three last terms are the corrections due to the cosmological constant. 
The classical advance in the perihelion is recuperated for zero cosmological 
constant (A — > 0) and the complete expression agrees with the one reported 
in |4|. This derivation is much simpler than the obtained in the work of G. 
Kraniotis and S. Whitehouse JS], where the evaluation was done by means of 
the genus 2 siegelsche modular forms and including the Mercury's data. 
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